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Abstract 

The detection of the Cosmic Microwave Background Radiation (CMB) was one of the most 
important cosmological discoveries of the last century. With the development of interferometric 
gravitational wave detectors, we may be in a position to detect the gravitational equivalent of 
the CMB in this century. The Cosmic Gravitational Background (CGB) is likely to be isotropic 
and stochastic, making it difficult to distinguish from instrument noise. The contribution from 
the CGB can be isolated by cross-correlating the signals from two or more independent detectors. 
Here we extend previous studies that considered the cross-correlation of two Michelson channels 
by calculating the optimal signal to noise ratio that can be achieved by combining the full set of 
interferometry variables that are available with a six link triangular interferometer. In contrast to 
the two channel case, we find that the relative orientation of a pair of coplanar detectors does not 
affect the signal to noise ratio. We apply our results to the detector design described in the Big 
Bang Observer (BBO) mission concept study and find that BBO could detect a background with 
Q gw > 2.2 x 10 -17 . 
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I. INTRODUCTION 


Most theories describing the formation of the Universe (Inflation, Qflation, etc.) predict 
that processes in the early Universe will lead to the copious production of gravitational 
waves. The detection of such a cosmic gravitational background (CGB) would allow us to 
probe the earliest moments in the history of the Universe, and place strong constraints on 


the competing theories [l|. However, detecting the CGB will not be easy since it will be 
hidden behind the signals from astrophysical sources (binaries systems...) and buried in the 
instrumental noise. Our current understanding of compact binary systems suggests that 
there is a window above 0.1 Hz where the number of astrophysical sources is small enough 
that their contribution can be isolated and removed from the detector data streams 


he instrument noise by cross-correlating the outputs 


The CGB signal can then be dug out of 

of two or more independent detectors Jj, fl], The Big Bang Observer (BBO) 6| has been 
proposed as a future space based mission designed to operate in the range 0.1 —^ 1 Hz. The 
BBO proposal calls for a fleet of triangular interferometers operating on the same principle 
as the Laser Interferometer Space Antenna (LISA). The BBO detectors will be ~ 100 times 
smaller than the LISA detector, and will be considerably more sensitive. It is possible 


to synthesize three independent data channels [7] (labeled ( A,E,T )) in each detector. In 
principle one could cross correlate these channels using data from a single detector as the 
channels are nominally noise orthogonal. However, the A and E channels are also signal 
orthogonal, and the T channel has poor sensitivity to waves with wavelengths larger than 
the detector. Moreover, the three channels are constructed from links that share common 
noise sources, so it will be difficult to achieve exact noise orthogonality in practice. The BBO 
design overcomes these obstacles by employing multiple detector units. A pair of co-planar 
detectors yields the greatest sensitivity as the antenna patterns have significant overlap. 
The BBO design also calls for two widely separated “outrigger” constellations that give 
enhanced angular resolution for detecting astrophysical sources fl , but the wide separation 
(v^AU = 866 sec) renders them useless for performing cross correlated detection of the 
CGB in the 0.1 —> 1 Hz range. 

Here we study the optimal cross correlation of co-planar triangular interferometers. This 
work generalizes earlier studies fly, q that only considered the cross correlation of two 


Michelson channels. By cross-correlating all possible combinations of the A, E, T channels 


2 






in the two detectors the overall sensitivity is improved by a factor of \[2 at low frequencies 


i quen 

B 


and by a factor of \/3 at high frequencies. In contrast to the single channel case jll|, the 
optimal sensitivity is independent of the relative angle between the two co-planar detectors. 
We find that the fiducial BBO design |]] will permit the detection of a scale invariant CGB 
with Q gw = 2.2 x 10~ 17 at 95% confidence. 

We begin in Section II by describing the gravitational wave response of the A, E and T 
channels. This is followed in Section III by a calculation of the noise transfer functions in 
each channel. Section IV describes the optimal cross correlation of the independent data 
channels in a pair of co-planar detectors. Section V covers the numerical evaluation of the 
overlap functions and a calculation of the optimal BBO sensitivity. We use geometric units 
with G — c — 1. 

II. DETECTOR RESPONSE 


The noise orthogonal data channels A, E and T are formed from linear combinations of 
the three Sagnac channels si,s 2 ,s 3 : 

A = -y^( s 3 ~ s i) E — -^=(si — 2 s 2 + S 3 ) T — -^=(si + s 2 + S 3 ) (1) 

The Sagnac interferometer measures the phase difference of two laser beams starting from 
the same location and going around the triangle formed by the three spacecraft, one trav¬ 
eling clockwise, the other counterclockwise. Ideally the phase difference is due only to the 
variations of the interferometer arms’ length caused by the gravitational waves. Therefore 
if the beams start from spacecraft 1 , the signal is simply 

Si(t) = t^-[Zi 3 (Z — 3L) + I 32 (t — 2 L) + Z 2 i(t — L) — Z 12 (t — 3 L) — Z 23 (t — 2 L) — Z 3 i(f — L)], (2) 

where (t — nL) is the distance at time t — nL between spacecraft i and j and L is the 
length of the interferometer arms (assuming all the arms have the same length). For a plane 
gravitational wave propagating in £2 direction, this can be shown [llj to reduce to 


si(t) = D s (ai, & 2 , c 3 , 0 , /) : h(/, t, .xq), 


( 3 ) 


where a\, & 2 , c 3 are vectors that point along the interferometer arms, h(/, t,x 1 ) is the tensor 
describing the wave in the transverse-traceless gauge at point x[, 


D(ai, & 2 , c 3 , £2, /) = - (a® a Ti(/,a) + b® b T 2 (f, b) + c® cT 3 (/,c)) 


( 4 ) 
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and 


Ti(a, f ) = e -^"( 1+ “^)sinc (/„( 1 + a • Cl)) - e^^sinc (/„(1 - a • ft)) , (5) 

T 2 (6, /) = e ~ i/n [ 3+(a-5 H sine (/ n ( 1 + b ■ ft)) - sine (/ n (l - b ■ ft)) , (6) 

T 3 (c, /) = e-^( 5 - c ^sinc (/ n (l + c ■ ft)) - e^-^sinc (/„( 1 - c ■ ft)) , (7) 

with f n = nLf. Similarly the Sagnac signal extracted at vertex 2 and 3 can be found from 
symmetry by rotating the system: 

s 2 (t) = D s (ci, a 2 , b 3 , ft, /) : h (f,t,x 2 ), 

s 3 (t) = D s (6i, c 2 , a 3 , ft, /) : h(/, t, x 3 ). (8) 

We are now ready to find the A, E and T detector responses to a plane gravitational 
wave. To simplify matters we write the signal in the form 

s n (t) = D n (ft, /) : h(/, t, x 0 ), (9) 

with 

D n (ft, f) = a®aT“ + b®bT b n + c®c T c n . (10) 

We define our reference point x$ to be the center of the triangle, and write 

h (f,t,x) = (11) 

The T variable, which is also called the symmetrized Sagnac, has an obvious cyclic symmetry: 

= Tf = T t , (12) 

where 

- e -i^(9+(wi-tr 3 )-o) , \ r 

T t {u 2 ■ ft, f) = --- [l + 2 cos(2 f n )J sinc(/ n (l + u 2 • ft)) 

-sinc(/ n (l - u 2 ■ ft)) (13) 

The variable A, on the other hand, does not have this nice symmetry: 

T%(a-Cl,f) = ^=sin(/„)e"^( 6+(< ^ 6) '^ sinc(/ n (l + a ■ Cl)) 

+e _2l/n sinc(/ n (l - a ■ ft)) , (14) 


4 



T\(b • ft, /) = «) [ e - 2 ^~.sinc(/„(l + 6 • ft)) 

+sinc(/„(l - 6 - ft)) 


( 15 ) 


I5(c- ft,/) = g-|sm(/ n ) e -‘^( 9+ < 5 - s > n ) 


sine (/n(l + C ' 
+sinc(/ re (l - c • 0)) 


(16) 


nor does the variable E: 


Tfi(a ■ Cl, f) — 


6VQL 


sme 


(/ n ( 1 + a • 0)) ( e -^( 9 +P-?)-n) + e -i£(3+M)-n) 


_2 e -^( i5 +(c-6).h)) + sinc (/ n (i _ s. n)) (^e-^^+iS-byn) 

_ e -i^(l5+(c-b) Cl) _ e -i^(9+(c-6)-n)^ 


(17) 


T b E (b-nj) = 


6^/6 


sine 


(/ n (l + b • ft)) ( e -^(i5+(s-5)-h) + e -<^(9+(ff-5)-n) 
-2 e -^( 3+{s - 5) ^)) + sinc(/„(l - 6 • ft)) (2e" i ^( 15+(s “ 5) ' A ) 

_ e -i I ^(9+(a-c)-Cl) _ e -iE.{ 3+(a-c)-S7) 


(18) 


T c E (c-n,f) = 


6^6 


sine (/ n (l + c- ft)) (e-^( 3+ ^ s ^) + e -^( 1 s+(M-n) 

_ 2e -^( 9 +(6-5)'h)) + ^(^(1 _ c • ft)) ^2e-^( 9+(g - s) '^) 

_ e -i^( 3 +(6-a)-o) _ e -i^(l5+(6-a)-n)^ 


(19) 


III. NOISE SPECTRAL DENSITY 

Until now we have only considered the gravitational wave contribution, to the 

time-varying part of the phase (j>ij(t). Our next task is to account for the instrument noise 
contributions. There are three main noise sources: the laser phase noise (7(f), the position 
noise n p {t) and the acceleration noise n a (t). 
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The total phase variation is given by 


= Ci(t - Lij) - Cj(t ) + -0 a(t ) + nfj(t) - x tj • [n“ (t) - n“ (t - L^)]. (20) 


The position noise (t) includes shot noise and pointing jitter in the measurement of 
the signal sent by spacecraft i and measured by the photo-detector in spacecraft j. The 
acceleration noise n“ (f) is from the gravitation reference system in spacecraft j along the axis 
that points toward spacecraft i. The phase noise associated with the laser on spacecraft i is 
denoted Ci(t). ft is easy to show that the phase noise cancels in a rigid, non-rotating Sagnac 
interferometer. More complicated second generation Sagnac variables can be constructed to 
account for the rotation and flexing of the array 1121 . Il3| . For simplicity we work with the 
basic Sagnac variables as they give results that are almost identical to those found using the 
second generation variables. 

We assume that all the interferometer arms are of approximately equal length (L = 
5 x 10 7 m), and that the noise spectral densities S n (f ) are similar on each spacecraft. The 
noise transfer functions are then given by 


sj(/) = 2 [1 + 2cos(2/„)f SJ(/) + 4sin 2 (/„)S“(/) 


( 21 ) 


where according to the fiducial BBO design 


9 0 x 10 -34 

S a n (f) - - Hz 


-l 


( 22 ) 


(2 I r/)-‘(3L) 2 

This expression for the Symmeterized Sagnac noise transfer function was previously derived 
in Refs. [7, [ll|. The noise transfer functions in the A and E channel share the same form, 
as first pointed out in Ref. [7]: 


Sn(f) = S n(f) = 8sin 2 (/ n ) 


(2 + cos(2 f n ))S„{f) + 2(3 + 2cos(2/„) + cos(4/ n ))S , “(/) 


(23) 


IV. CROSS-CORRELATION OF TWO DETECTORS 

The CGB signal can be extracted from the instrument noise by cross-correlating the 
outputs of two independent interferometers. The pair of co-planar interferometers of the 
BBO do not share any common components, so the noise in each detector should be largely 
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uncorrelated. Possible correlated sources of noise include solar flares and fluctuations in the 
refractive index of the inter-planetary medium. Another potential source of correlated noise 
is the residual from subtracting foreground sources such as double neutron star binaries US]. 
Here we will assume that any correlated sources of noise are well below the level of the CGB. 

We assume that the CGB is stationary, Gaussian, isotropic, and unpolarized. The back¬ 
ground can be expanded in terms of plane waves: 

M*,2) = E r d f (24) 

I J — oo j 

where e" 4 (H) are polarization tensors given in term of the basis tensors 

e + (0, VO = e + (fi) cos 2 if; — e x (fi) sin 2^, 

e x (f2,V0 = e + (fi) sin 2^ + e x (fi) cos2-0. (25) 

The basis tensors can be expressed in term of an orthonormal set of unit vectors rh, n and 
hi as 

e + (fi) = 

e x (H) = in ® fi + n ® rh. (26) 

From the strain 

S(Cl, /, t) = D (Cl, f ) : ( h + (f , t, x 0 )e+(n, 4>) + h* (/, t, X 0 )e x (Cl, ^)), (27) 

we hnd that after averaging over polarizations, the cross-correlated signal of two detectors 
in the presence of such a background is given by 

/•OO 

(Si(t)s 2 (t)) = / dfS h (f)R 12 (f), (28) 

Jo 

where Sh(f) is the total power spectral density due to both polarizations: 

S k {f) = St(f) + SZ(f), (29) 

and 

RvU) = E / A /)■ (30) 

The transfer function R\ 2 is a purely geometric factor that accounts for the overlap of the 
antenna patterns of the two detectors. The antenna pattern functions F are given by 

F A (Cl,f) = B(n,f):e A (Cl). (31) 
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The optimal signal-to-noise ratio squared is given by 

\R 12 (f)\ 2 


SNRp = 2 T 


dfS z h (f)- 


(32) 


S ni (f)Sn 2 (fY 

where S ni is the noise spectral density of interferometer 1, S n2 is the noise spectral density 
of for interferometer 2 and Sh is the spectral density of the CGB. A signal to noise ratio of 
SNR = 3.3 indicates that the CGB has been detected at 95% confidence, with a 5% false 
alarm probably Q. 

The power spectral density of the CGB is related to the energy density in gravitational 
waves per logarithmic frequency interval, Q gw (f) (in units of the critical denisty), by 




(33) 


47T 2 / 3 

Standard inflationary models predict that Q gw (f) will be roughly scale invariant, with an 
amplitude Q gw ~ ICC 15 in the / = 1 Hz region. 

The contribution to the cross-correlated SNR per logarithmic frequency interval can be 
written as 

dSNR 2 c K pt (f ) 

din/ K s (f) 1 J 

where h e s(f) is the effective sensitivity curve 


M/) = 


Sn(f) 


\ |*12(/)| 

and h op t(f) is the optimally filtered CGB signal 

Ao P «(/) = (2T/) 1 ' 4 v /%(7) 


(35) 


(36) 


These definitions are equivalent to the usual definitions used to plot sensitivity curves and 
gravitational wave signals for coherent gravitational wave sources. The main difference is 
that the optimally filtered signal strength grows as T 1 / 2 with coherent matched filtering, 
while it only grows as T 1 / 4 for cross-correlated stochastic signals. 


V. NUMERICAL ANALYSIS 

The optimal SNR for the BBO comes from combining the full set of independent inter¬ 
ferometry channels in each of the two co-located detectors that form the star constellation: 

snr4, = £ E SNR C ( 37 > 

ot=A\,E\,T\ (3=A 2 ,E2,T2 
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Our first task is to calculated the geometrical overlap of each pair of channels, R a g[f). 
The all-sky integral in was performed numerically using the HEALPIX package jl5| . 
Plots of the overlap factors are shown in Figures 1-6 for the standard BBO configuration 
in which the two overlapping detectors form a symmetric six pointed star. The plots show 
the R 1 s scaled by overall factors such as sin 2 (/ n ), which they share in common with the 
corresponding noise transfer functions. 



We expect R 12 to vanish for the cross terms as the A, E and T channels are approximately 


signal orthogonal. I 11 1 
to Cutler’s sj and sjj 


to be equivalent 


re low frequency limit the A and E can be shown 

variables, which describe two Michelson detectors rotated by an 
angle 7 t/ 4. Combining the geometrical factors R 12 with the noise transfer functions leads to 
the effective sensitivity curves shown in Figures 7 and 8. 

The combined effective sensitivity follows from the optimal signal to noise ratio (ED • 
Figure 9 shows the combined sensitivity curve using all channel combinations plotted against 
the optimal CGB signal for a scale invariant spectrum with Q gw = 10~ 15 . 

An alternative way of conveying the information contained in Figure 9 is to plot SNR(/) = 
(dSNR 2 (/)/dln/) 1 / 2 : the contribution to the signal to noise ratio per logarithmic frequency 
interval. The optimal SNR(/) for Vt gw = 10~ 15 is shown in Figure 10. 

The preceding graphs assumed that the two triangular interferometers formed a symmet¬ 
rical six pointed star. However, one could consider alternative arrangements where the two 
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detectors are rotated by an arbitrary angle A with respect to one another. In our numbering 
convention for the vertices of each triangle the symmetric star corresponds to a rotation 
angle of A = n. As we vary A, the individual SNRs from A { x A 2 and E\ x E 2 decrease, 
with minima at A = 7 t/ 4 and A = 37 t/ 4 as expected 11|. The SNR for Ai x E 2 , on the 
other hand, increases, with maxima at A = 7 t/ 4 and A = 37 t/ 4. This behavior is illustrated 
in Figures 11 and 12. The net effect is that the optimal SNR is independent of the relative 
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FIG. 4: R\ 2 {fn)/ ((1 + 2cos(2/ n ))sin(/ n )) for A 1 x T 2 or T, x A 2 



orientation A. 

We find that the optimal SNR for scale invariant CGB that can be achieved by the fiducial 
BBO design is equal to 


SNR opt = 155 1 


t n 


gw 


Hn 


5yr 10 15 \ 70 kms 1 Mpc 


-i 


(38) 


which is well above the 3.3 threshold mentioned in the introduction. Conversely, the 
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FIG. 6: Ruifn)/ sin 2 (/ n ) for A ± x E 2 or E 1 x A 2 . 



FIG. 7: Effective sensitivity for A\ x A 2 ( E\ x E 2 ) 


minimum Q gw for which we could detect the CGB with the 95% confidence is 


2.2x10 17 (5yr/T) 1//2 (70kms 1 Mpc 1 /H 0 ) 2 . Our findings agree with the recent 
dent calculation by Seto M- 


equal to 
indepen- 
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FIG. 8: Effective sensitivity for T\ x T 2 
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FIG. 9: The combined sensitivity curve (solid line) plotted against the optimal CGB signal (dashed 
line) for an observation time of T = lyr. 

VI. CONCLUSION 

We have determined that the fiducial BBO design would be able to detect a scale invariant 
CGB with an energy density as low as = 2.2 x 1CT 1 ' (5yr/T) 1 / 2 (70kms _1 Mpc _1 /-^o) 2 - 
We found that the optimal sensitivity is independent of the relative orientation of the co- 
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FIG. 10: Optimal SNR(/) 



FIG. 11: SNR Al A 2 (f) for an angle of A = n (dashed line) and A = 7r/4 (solid line), 
planar detectors used to perform the cross correlation. 
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